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Abstract. For commuting linear operators Pq, Pi, ■ ■ ■ , Pi we de- 
scribe a range of conditions which are weaker than invertibility. 
When any of these conditions hold we may study the composition 
P = PqPi ■ ■ ■ Pi in terms of the component operators or combi- 
nations thereof. In particular the general inhomogeneous problem 
Pu = f reduces to a system of simpler problems. These prob- 
lems capture the structure of the solution and range spaces and, if 
the operators involved are differential, then this gives an effective 
way of lowering the differential order of the problem to be studied. 
Suitable systems of operators may be treated analogously. For a 
class of decompositions the higher symmetries of a composition P 
may be derived from generalised symmmetries of the component 
operators Pi in the system. 



1. Introduction 

Given a vector space V and a system P , Pi, • • • ,Pe of mutually com- 
muting endomorphims of V we study the composition P := P$P\ ■■■ Pe- 
lt is natural to ask whether we can reduce the questions of null space 
and range of P to the similar questions for the component operators 
Pi. If these component operators are each invertible then of course 
one trivially has a positive answer to this question. On the other hand 
experience with, for example, constant coefficient linear ordinary differ- 
ential equations shows that this is too much to hope for in general. Here 
we review, discuss, and extend a recent work [6] in which we introduce 
a range of conditions which are significantly weaker than invertibility 
of P and yet which, in each case, enables progress along these lines. 

Each condition we describe on the system P , Pi, • - • , Pi is termed an 
a-decomposition (where a is a subset of the power set of the index set 
{0, ■ • • ,£}). The case that the operators Pj are each invertible is one 
extreme. Of course one may ask that some of Pj are invertible but, ex- 
cluding an explicit assumption along these lines, the next level is what 
we term as simply a decomposition. This is described explicitly in Sec- 
tion [2] below, but intuitively the main point is that each pair Pj and 
Pj, for i ^ j in {0, • • • ,£}, consistes of operators which are relatively 
invertible in the sense that for example Pj is invertible on the null space 
of Pj and vice versa. In the case that we have a decomposition then 
one obtains very strong results: the null space of P is exactly the direct 
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sum of the null spaces for the factors P,; the range of P is precisely the 
intersection of the range spaces for the factors; and one may explicitly 
decompose the general inhomogeneous problem Pu = f into an equiv- 
alent system of "lower order" problems Pj«j = /, i = 0, • • • , £. (Note 
that the fact the same inhomogeneous term / 6 V appears in Pu = f 
and in each of the pWj = / problems is one signal that the construction 
we discuss is not the trivial manouever of renaming variables.) 

At the other extreme of the a-decompositions we ask only that the 
operator (P , • ■ ■ , Pe) : V — ► ® l+1 V is injective with left inverse given 
by a system of V endomorphisms Qi, i = 0, ••• ,£, which commute 
with the PjS. Remarkably this is sufficient for obtaining results along 
a similar line to the case of a decomposition, but the extent of sim- 
plification is less drammatic: issues of null space and range for P are 
subordinated only to similar questions for the operators P % := P/p. 
The full summary result for a-decompositions is given in Theorem 12.11 
below. 

A natural setting for the use of these results is in the study of opera- 
tors P which are polynomial in a mutually commuting system of linear 
operators T>j : V — > V, j = 1, • • • , k. This is the subject of Section 
[HI Given a P which is suitably factored, or alternatively working over 
an algebraically closed field, one sees that generically some algebraic 
a-decomposition is available. The main point here is the word "alge- 
braic". The a-decompositions of compositions P = PoPl ■ ■ ■ Pe involve 
identities which involve operators Qi, i — 0, • • • ,£, which invert some 
subsystem of (Po, • ■ ■ ,Pe)- In the case of an algebraic decomposition 
the QiS are also polynomial in the VjS. So, for example, if the factors 
Pj of P are differential and polynomial in the differential operators 
T>i then the Qi are also differential and are given in terms of the T>i 
by explicit algebraic formulae. In particular pseudo-differential calcu- 
lus is avoided. These results are universal in the sense that they are 
independent of any details of the operators Vj. 

The idea behind the decompositions of the equation P • • • P(U = f is 
rather universal; for example one can extend it to systems of equations. 
In this direction our aim is mainly to demonstrate the technique, so we 
shall treat, as an example (see Section 0J), only one specific situation 
where the idea applies. This is a system of k + 1 equations where the 
first equation Pq - ■ ■ Pgu = f is factored and the remaining ones are of 
the form R^u = g\ 1 < j < k for given /, gi G V and G End(V). 
This may be viewed as a problem where one wants to solve the problem 
P • • ■ P(U = f, subject to the conditions Ry'u = gi . The difference, in 
comparison to the single equation problem, is that now the operators 
i?w feature in the relative invertibility of the factors Pj of P. The 
result is that provided one has a suitable decomposition at hand, the 
original system is equivalent to a family of "lower order" systems of the 
same type as the original one. 
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Finally in Section [5] we discuss symmmetries of operators. We define 
a formal symmetry of an operator P : V — > V to be an operator 
S : V — > V such that PS = S'P for some other operator S' on V. For P 
a Laplacian (or Laplacian power) type operator differential operator, 
and S, S' differential, such symmetries are central in the separation 
of variable techniques [H IS]- For P = PqP\ ■ ■ ■ Pi, as above, the tools 
we develop earlier are used to show that the formal symmetry algebra 
of P is generated by the formal symmetry operators, and appropriate 
generalisations thereof, for the component operators PqP\, ■ ■ ■ ,Pg. 

The first author would like to thank the Royal Society of New Zealand 
for support via Marsden Grant no. 06-UOA-029. The second author 
was supported from the Basic Research Center no. LC505 (Eduard 
Cech Center for Algebra and Geometry) of Ministry of Education, 
Youth and Sport of Czech Republic. 

2. Decompositions and ^-decompositions 

Let V denote a vector space over a field F and consider linear oper- 
ators Pq, . . . , Pi (i.e. endomorphisms of V) which mutually commute. 
In [6] we study properties of the operator 

(1) P = P P 1 ---P e . 

For example, an obvious question is: what can we say about the ker- 
nel and the image of P in terms of related data for the component 
operators P{t This is clearly straightforward if the operators Pi are 
invertible, but the point of our studies is that much weaker assump- 
tions are sufficient to obtain quite striking results. These assumptions 
are captured in the notion of various "decompositions"; the different 
possible decompositions are parametrised by a nonempty system of 
subsets of L := {0, . . . ,£}. We shall use the notation Pj := YljejPj 
for 7^ J C L and we set Pq := idy. Further we put P J := Pl\j, 
Pi := P^} and write 2 L to denote the power set of L. Also \ J\ will 
denote the cardinality of a set J. 

Definition. For a linear operator P : V — > V, an expression of the 
form ([]]) will be said to be an a -decomposition of P, with 0^«C 2 L , 
L ^ a, if there exist operators Qj e End(V), J G a such that 

(2) id v = ®J pJ ' [Pi,Pk] = [Qj,Pi\ = 0, i,keL,Jea. 

The choice a := {{0}} means that P (hence also each of the Pj) 
is invertible. The other possible decompositions involve weaker as- 
sumptions on the component operators. At the next level is the en- 
decomposition with a := {J C L \ \J\ = 1} which will be termed 
simply a decomposition of P. In this case we still obtain, for example, 
that M{P) = 0A/"(Pj). Therefore the problem Pu = 0, u G V is re- 
duced to the system = 0, Ui G V for i G L. In the case that the p 
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are differential operators, this result shows that, given a decomposition, 
the equation Pu = reduces to the lower order system p-Uj = 0. 

For the general a-decomposition we do not generally obtain a direct 
sum analogous to N(P) = ®A/"(Pj) as above, however we still get a 
reduction to a "lower order" problem. The key is the following theorem 
which is a central result in [6]. (See the latter for the proof and more 
details.) 

Theorem 2.1. [6] Assume P : V — > V as in is an a-decomposition. 
Let us fix f G V. There is a surjective mapping B from the space of 
solutions (uj)j£ a G ®' Q 'V of the problem 

(3) Pjuj = /, Jea. 

onto the space of solutions u G V of Pu = f . 

Writing Vp for the solution space of Pu = f and (for J G a) Vj for 
the solution space of PjU = f. The map B : Xj ea Vj — > Vp is given 
by 

JGo 

t4 right inverse for thts is F : V f P -> X JgQ V^ <?iwen (component-wise) 
by 

u i— > P j m ; 

on V we /iawe B o F = idy. 

If a satisfies I H J = 0, /or all I ^ J E a, then F is a 1-1 mapping 
and F o B is the identity on the solution space to 

Remark 2.2. 1. The important feature of the decomposition of inho- 
mogeneous problems is that in Pu = f and (|3j) it is the same / G V 
involved. So (J3j) describes the range P, 1Z(P), in terms of the range of 

the p^ n(P) = ni =Q n{Pi). 

2. The condition in the last paragraph in the theorem is satisfied by 
a decomposition, but is easy to construct other examples. In any case 
where this is satisfied the mappings F and B are bijections and we do 
get a direct sum decomposition of M{P). The point, which is easily 
verified, is that from j2j) it follows that for each Jea 

QjP J : Af(P) - Af(Pj) 

is a projection. 

Although the a-decomposition (J2J) is what is directly employed in the 
previous theorem and its proof, there is a distinct, but related, notion 
which shows what it really means for the commuting operators Pj. The 
following definition introduces an idea of a decomposition which turns 
about to be in a suitable sense "dual" to the previous one. 

Definition. For a linear operator P : V — > V, an expression of the 
form (pQ) will be said to be a dual (3 -decomposition of P, ^ (3 C 2 L , 
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{0} ^ (3 if for every J G (3 there exist operators Qjj G End(V), j G J 
such that 

(4) id v = Y,Qjj p j> [Pi,Pk} = [Qjj,Pi\ = 0, i,keL,jeJ. 

To describe the suggested duality (see Proposition 12.41 below), first 
observe that each system a C 2 L is partially ordered by restricting 
the poset structure of 2 L . The sets of minimal and maximal elements 
in a will be denoted by Min(a) and Max(a), respectively. We say the 
system (3 C 2 L is a /ower setf, if it is closed under taking a subset. (That 
is, if I G /3 and J C I then J G The upper set is defined dually. The 
lower set and upper set generated by a system a C 2 L will be denoted 
by £(a) := { J C / | / G a} and W(a) := { J D / | J C L and / G a}, 
respectively. 

The proof of the following is obvious. 

Lemma 2.3. Let a C 2 L . T/ien P = P - ■ ■ Pg satisfies the following: 

(i) it is an a- decomposition •<=>- is a Max(a) -decomposition -^=> it 
is an C(a) -decomposition 

(ii) it is a dual a -decomposition <^=^ it is a dual Min(a) -decomposition 
<^=^> it is a dual U (a) -decomposition. 

To formulate the relation between a- and dual a-decompositions, we 
need the following notation. We put a u := 2 L \C(a) and a 1 := 2 L \U(a). 
Clearly (a u ) 1 = C(a) and (a l ) u = U(a). Also it is easily seen that 

a u = {J C L | V/ G a : J\I ^ 0} 
a ! = {JCI|V/Ga:/\J^0}. 



(5) 



Proposition 2.4 (The duality). (G|) is an a -decomposition if and 
only if it is a dual a u -decomposition. Equivalently, (Qp is a dual (3- 
decomposition if and only if it is a (3 l -decomposition. 

In particular, is a decomposition if and only if it is a dual (3- 
decomposition for (3 := {J C L \ \ J\ — 2}. This means 

(6) id v = QijPi + Qj,iPj 

where Qi.j G End(V) and satisfy [Qij, Pk] = for every triple of inte- 
gers k) such that < i, j, k < £ and i j. 

When (P) is satisfied we shall say that the operators Pj and Pj are 
relatively invertible. The dual version of a (true) decomposition is the 
dual /^-decomposition for (3 = {J C L \ \ J\ = 2}; this will be termed 
simply a dual decomposition. The general dual /^-decomposition means 
that for every J G (3, the operators Pj, j G J are relatively invertible. 

Remark 2.5. 1. The proof of Proposition 12.41 in [6] is constructive in 
the sense that starting with an a-decomposition (j2J), there is a simple 
recipe which describes how to construct the operators Qjj from (pj , as 
required for the dual (^"-decomposition, and then vice versa. 
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2. The operators Qj in (J2j) and Qjj in (jlj) are not given uniquely. 
Also we can see an obvious duality between Qij, Qjj and Pj, Pj in 
(jHJ). This is discussed in [6 J where notion of (dual) a-decompositions 
is formulated in the language of Koszul complexes. 

From the practical point of view, given an operator P : V — > V as 
in (JTJ), to apply Theorem 12.11 one needs to show whether P is a (dual) 
a-decomposition and also to determine explicitly the corresponding 
operators Qj (or Qjj in the dual case). Also, one can ask which choice 
of a yields the most suitable a-decomposition. Another strategy might 
be to "regroup" the operators Pi (e.g. to consider the product pPj as 
a single factor) and then to seek a better a-decomposition. In the 
case that the operator P is polynomial in other mutually commuting 
operators T> } ■ : V — > V, there is a category of decompositions which 
arise algebraically from the formula for P. Within this category all 
these questions can all be solved in a completely algorithmic way. 

3. Operators polynomial in commuting endomorphisms 

and algebraic decompositions 

Writing V to denote a vector space over some field F, suppose that 
T>i : V — > V, i — 1, • • • , k, are non-trivial linear endomorphisms that are 
mutually commuting: T>{Dj = T>jT>i for i, j G {1, •••,&}. We obtain a 
commutative algebra ¥[T>] consisting of those endomorphisms V — > V 
which may be given by expressions polynomial (with coefficients in F) in 
the T>i. We write x — (x%, . . . , %k) for the multivariable indeterminate, 
and ¥[x] for the algebra of polynomials in the variables x%, . . . , x^ over 
the field F. There is a unital algebra epimorphism from ¥[x] onto ¥[T>] 
given by formally replacing each variable Xi, in a polynomial, with T>j. 

The simplest case is when k — 1, that is operators polynomial in a 
single operator V. We write ¥[D] for the algebra of these. Since any 
linear operator T> : V —>■ V is trivially self-commuting there is no re- 
striction on T>. Thus this case is an important specialisation with many 
applications. In this setting we may quickly find algebraic decomposi- 
tions. Let us write ¥[x] for polynomials in the single indeterminate x 
and illustrate the idea with a very simple case. Consider a polynomial 
P[x] = (x + A ) (x + Ai) • • • (x + Xi) where for i = 0, • • • , £, the Xi e ¥ 
are are mutually distinct (i.e. i ^ j : =^> Aj ^ Xj). Related to P[x] are 
the polynomials obtained by omitting a factor 

i 

p>}= n ( x+x ^- 

Then we associate to P[x] the following decomposition of the unit in 
¥[x}. 
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Lemma 3.1. 

1 = a P°[x] + aiP 1 ^] + • • • + a t P l \x\ , 

where 



Aj — A,; 



1 



i^j=0 



Proof. For £ = we take the first display to mean 1 = 1. For £ = 1 the 
Lemma states that 

1 = r 1 r (x + Ai) + 1 (x + A ) 

Ai — Ao Ao — Ai 

which is clearly true. Now assume that the result holds for all polyno- 
mials with £—1 factors. In particular associated to P l [x] = Y[iZo( x +\i) 
and P°[x] = Yli=i( x + Ai) we have the identities 

(7) 1 = /3 Q [x]-\ h/ViQf-iM and l = 7 1 i2 1 [x]H \-j e Ri[x], 

respectively, where we have 



Ao — Ai 



and 



£ 1 \ 

Rk[x]= JJ (x + Aj), 7fc = II T — ^ , for fc = l, 

Now multiplying the first identity of by (x + A^) , multiplying the 
the second identity of by (x + A ) and then taking the difference 
yields 

A, - Ao = f3 P [x] + (Pi - ll )P 1 [x] + - ■ ■ + (&_x - li . 1 )P t - 1 [x)-'y e P e [x). 

This establishes the result as Po/(^£ — Ao) = cuo, "fe/(\e — Ao) = — ae, 
while for i — 1, • • • , £ — 1 we have 

A ~7i = / TT 1 w 1 1 ^ 1 = ^ 

A^ — Ao 7" , A 7 - — Aj Ao — Aj A^ — A« A^ — Ao 

□ 

Thus we have the following. 

Proposition 3.2. For P[D) = (V + A )(£> + Ai) • • • (V + X e ) we have 
a decomposition given by 

td v = Q P°[V} + --- + Q f P e [V}, 
where P l = H ijij=0 (V + Xj) and Q, t = ai for i = 0, • • • ,£. 
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Thus we may immediately apply Theorem l2.ll and in fact the stronger 
variants for decompositions as in |6j, to reduce homogeneous or in- 
homogeneous problems for R to corresponding problems of the form 
(V + X i )u i = f. 

Remark 3.3. However the point we wish to emphasise heavily is that 
we used no information about the operator V to obtain the decomposi- 
tion in Proposition [321 D can be any linear operator V : V — » V on any 
vector space V. Thus we will say that Proposition 13.21 is an algebraic 
decomposition of -P[X*]. For specific operators T> there may be other 
decompositions (or a-decompositions) that do use information about 
V. 

For operators polynomial in mutually commuting operators T>{ we 
generically may obtain a-decompositions that are algebraic in this way; 
that is they arise, via the algebra epimorphism ¥[x] — > ¥[D], from a 
polynomial decomposition of the unit in ¥[x]. These are universal 
a-decompositions that are independent of the details of the T>i, i — 

l,--- ,k. 

Via the Euclidean algorithm, and related tools more powerful for 
these purposes, we may easily generalise Lemma \3. II to obtain decom- 
positions for operators more interesting than P[D] as in the Proposition 
above. The case of a operators polynomial in a single other operator is 
treated in some detail in [6] so let us now turn our attention to some 
general features which appear more in the multivariable case k > 2. 

Given polynomials -Pot*], AWr" > Pt\ x ] £ ¥[x] consider the prod- 
uct polynomial 

(8) P[x] = P [x]P 1 [x]---P e [x]. 

With L — {0, 1, • • • , L}, we carry over, in an obvious way, the labelling 
from Section [2] via elements of the power set 2 L ; products of the polyno- 
mial Pi[x] are labelled by the corresponding subset of L. For example 
for JCL, Pj[x] means rije./ ^J'M' while P J [x] means P L \j[x]. 

Considering the dual /3-decompositions, we need to verify that for 
each I G (3 we have 

(9) 1 G (Pi[x] : i G /) 

where (..) denotes the ideal in ¥[x] generated by the enclosed polyno- 
mials. It is useful to employ algebraic geometry to shed light on this 
problem, in particular to use the "algebra - geometry dictionary", see for 
example % Chapter 4]. Let us write M{S[x\) := {x G ¥ k \ S[x] = 0} 
for the algebraic variety determined by the polynomial S[x}. The ideal 
(Pi[x\) corresponds to the variety J\f% := M(Pi[x]) and the previous 
display clearly requires f] i€l Afi = 0. In fact if F is algebraically closed 
then the latter condition is equivalent to (J9j) . (This follows from the 
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Hilbert's Nullstellensatz, see [2].) Since generically f] ieI J^fi has codi- 
mension |/|, we conclude that (for F algebraically closed) if 1 1\ > k+1 
then in the generic case (J9]) will be satisfied. 

(Dual) decompositions and a- decompositions 

Aside from invertible P, the decompositions are the "best possible" 
among all a-decompositions (and similarly for the dual versions). How- 
ever they require 2 > k + 1 in the generic case (we need = 2 in (|§J)) 
which holds only for one variable polynomials. On the other hand 
there is always a chance that we obtain a decomposition by a suitable 
"regrouping" of the polynomials in ((HI) - So we can proceed as follows. 

Any polynomial P[x] G ¥[x] can be decomposed into irreducibles. If 
we were to take Pi[x] in ((8]) as such irreducibles then j8]) would not be 
generally the decomposition in the multivariable case. To obtain the 
decomposition one can consider products Pi[x] as single factors in ([8]) 
for suitable KL. This reduces the number of factors (i.e. £); to find 
an optimal (i.e. with £ maximal) version of this procedure we use the 
following lemma. 

Lemma 3.4. (i) Assume P[x\ has the form ([HP satisfying 1 G (Pi [a;], Pj[x 
for all < i < j < £ and P[x\ = Ro[x] ■ ■ ■ R r [x] is a decomposition 
of P[x] into irreducible polynomials Ri[x). 1/1 ^ (R p [x] , R q [x]} for 
some < p, q < r then there exists < i < £ such that Pi[x] = 
Rp[x]R q [x]P-[x] for a polynomial P/[a;]. 

(ii) Assume the polynomials Sq[x], . . . , S s [x] andT [x], . . . ,T t [x] sat- 
isfy 1 G (Si[x],Tj[x]) for all < i < s and < j < t . Then 
1 G (S[x],T[x\) where S[x] = S [x] ■ ■ ■ S s [x] andT[x] = T [x] ■ ■ -T t [x]. 

Proof (i) Assume the case Pi[x] = Rp[x]P"[x] and Pj[x] = R q [x}P"[x] 
for some i ^ j. Then 1 G (Pi [a;], Pj[x]) implies 1 G (R p [x], R q [x]}. 

(ii) We use the induction with respect to s + 1. Clearly the lemma 
holds for s + t = so assume s + t > 1. Then e.g. t > 1 so by 
the inductive hypothesis we get 1 G (5 [a;], T[x]) and 1 G (5[a;], T t [a;]) 
where T[x] = T [x] ■ ■ -T t -i[x\. This means 

1 = a[x]S[x] + b[x)f[x) and 1 = c[x)S[x) + d[aj]r t [aj] 

for some polynomials a[x], b[x), c[x] and d[x]. Now multiplying the 
right hand sides of these two equalities and using T[x] = T[x]T t [x], 
the lemma follows. □ 

We will use this lemma as follows. We start with the decomposition 
of P[x] = Ro[x] ■ ■ ■ R r [x] into irreducibles. Consider the graph with 
vertices t>o, . . . , v r , and an edge {v p , v q } for every < p, q < r such that 
1 ^ (R p [x], R q [x\). Denote the number of connected components by 
£+1 and the set of vertices in the ith component by G i7 < i < £. We 
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put 

P t [x] := J] R u [x], i = 0,... ,£ 

which yields the form (JSj) of P[x\. This satisfies 1 G (Pi[x), Pj[x\) for 
all < i < j < £ according to Lemma 13.41 (ii) and thus j8]) is the de- 
composition. Moreover, it follows from the part (i) of the lemma that 
no form P[x] = Pq[x] ■ ■ • P[,\x\ with £' > £ can satisfy the condition 
1 G (Pi[x], Pj[x}) for all < i < j < £'. (The discussed graph has £+ 1 
connected components and we need to "regroup" the vertices (corre- 
sponding to irreducible components) into £' + 1 groups corresponding 
to £' + 1 polynomials Pl[x}. If £' > £ then there is a pair of irreducible 
polynomials R p [x],R q [x] such that 1 G" [R p [x], R q [x]) which satisfy 
that R p [x) is a factor of P([x] and R p [x] is a factor of Pj[x] for some 

< i < j < £' ■ This is a contradiction with Lemma l3~4l (i).) 

Remark 3.5. From the geometrical point of view, if (jHJ) is a decom- 
position then M := M{P) = |JA/"j is the disjoint union. The previous 
paragraph describes how to find such decomposition for the variety M 
corresponding to any P[x] given by ©. Moreover, the obtained de- 
composition is minimal in the sense that Mi in M = \JMi cannot be 
disjointly decomposed into smaller (nonzero) varieties. 

Generically, the (dual) decompositions are not available in the multi- 
variable case. The "optimal" choice among all possible a-decompositions 
(in the sense of [6]) is as follows. The subsets ft C 2 L are partially or- 
dered by inclusion (i.e. now we use the poset structure of 2 2L ). Given 
an operator P in the form ([Sj) consider the family T of systems ft such 
that (jHj) is a dual /^-decomposition. Then T has the greatest element 
ftp = {Jptz-pft- Then an "optimal" choice for the dual /^-decomposition 
of P is ft := Mm(ftp). (We want to have in ft to the smallest possible 
subsets of L. So if the P^s are not invertible then the case of a dual 
decomposition may be regarded as the best we can do. With this phi- 
losophy we thus take ftp. Then using Lemma [2731 we take ft := Min(ftp) 
as it is easier to work with a smaller number of subsets.) Consequently, 
we obtain the optimal choice a := Max((/3p)') for the a-decomposition 
of P. 

Algorithmic approach and the Grobner basis 

Summarising, starting with P[x], the problem of obtaining a fac- 
toring ((8]) which is the decomposition or a suitable a-decomposition 
boils down to testing the condition 1 G (Pi[x\ | i G I) for various sub- 
sets I C L. This can be done using Buchberger's algorithm which 
computes a canonical basis (for a given ordering of monomials) of 
the ideal (Pi[x\ \ i G /), a so called reduced Grobner basis |2J. If 

1 G (Pi[x] | i G I), this basis has to be {!}. 
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In practice for reasonable examples this algorithm may be imple- 
mented in, for example, Maple. Actually, one can save some computa- 
tion and moreover obtain the explicit form of the operators Qj from 
(J2J) or Qjj from (jlj) by using Buchberger's algorithm without seeking 
the reduced basis. Consider the ideal I := (G) C F[x] (for a set of 
polynomials G) such that 1 6 I, and a Grobner basis G' of I. Then 
a G G' for some scalar a G F. The Buchberger's algorithm starts with 
G and builds G' D G by adding various linear combinations (with co- 
efficients in F[as]) of elements from G. So when a is added (and the 
algorithm stops) , it has the required form which expresses 1 as a linear 
combination of elements of G (up to a scalar multiple a). 

Example 

We shall demonstrate the previous observations on the operator 
P[D U V 2 ] :=V\V 2 + V\V\ + 3VfV 2 + £> 4 + W\V 2 2 + 3VfV 2 
+ 2V\ + 3V\V\ + V\V 2 + V{D\ - Pi. 

This correspond (after factoring) to the polynomial 

(10) P[x, y] = (x + l)(xy + y + l)x(x 2 + xy + x + y — 1). 

For example, taking V 1 — and V 2 = ^- the differential operator 
P[£,-§jj] ■ C7°°(M 2 ) -> C°°(M 2 ) is a the sixth order differential oper- 
ator. We apply the previous observation and Theorem 12.11 to reduce 
the corresponding differential equation to a lower order problem. In 
general, we start with the equation P[Di, T> 2 ]u = f for a given / G V. 

First we shall find the optimal (in the sense as above) dual (3- 
decomposition, /3 C 2^ 0,1 ' 2,3 ^ and/or whether we can obtain the dual 
decomposition after an appropriate regrouping of the factors in ( ITOl ). 
Many steps can be done directly in Maple. The factors 

Po[x,y] = x+l,Pi[x,y] = xy+y+1, P 2 [x, y] =x,P 3 [x,y] = x 2 +xy+x+y-l 

are irreducible; this can be verified by the command IsPrime. Using 
gbasis we see that 

1G (P [x,y},P t [x,y}), ie {1,2,3}, 

(11) l^{P^x,ylP j [x,y\), i,je {1,2,3} and 
1G (P 1 [x,y],P 2 [x,y],P 3 [x,y]). 

From the first two lines in (fTTT) and using the observation around 
Lemma IBT41 we conclude that P[x,y] = Po[x,y]P[x,y] is the decom- 
position for P[x, y] = (xy + y + l)x(x 2 + xy + x + y — 1). Following the 
first line in (fTTl) . one easily computes 

1 = -yP [x,y] + Pi[x,y], 1 = -(x - l)P [x,y] +xP 2 [x,y], 
l = (x+ y)P [x, y] - P 3 [x, y) 
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and multiplying these three relations we obtain 

1 = Q[x,y]P Q [x,y] +Qo[x,y]Pi[x,y]P 2 [x,y]P 3 [x,y] 

V v ' 

P[x,y] 

together with explicit form of the projectors Q and Q . Passing to 
the corresponding operators on the space V, this is the decomposition 
of P[V 1 ,V 2 ] = P [V U V 2 ] P[V U V 2 ]. Now using Theorem O (and 
Remark 12.21) we see that every solution u G V of P[Vi,V 2 ]u = f can 
be uniquely expressed as 

u = Q[V 1 ,V 2 \u + Q [V 1 ,V 2 ]u 

where u and u satisfy P [Vi,V 2 ]u = f and P[Vi,V 2 ]u = f. So we 
have reduced the original problem P[D 1 ,V 2 ]u = f to the system of 
latter two equations. 

Using the last line in (fTTi ). we can apply Theorem 12.11 to the equa- 
tion P[Vi,V 2 ]u = f. It is easy to compute the corresponding dual (3- 
decomposition for the operator P[V U V 2 ] = P 1 [V l ,V 2 ]P 2 [V 1 ,V 2 ]P ? ,[V l ,V 2 }, 
(3 = {{1, 2, 3}}; on the polynomial level we obtain 

1 = -Pi[x,y] + -(x + l)P 2 [x,y] - -P 3 [x,y]. 

This is actually also the a-decomposition, a = {{1, 2}, {2, 3}, {1, 3}} 
as P 1 [V 1 ,V 2 ] = P^[V X ,V 2 ] etc. according to the notation in 
Now applying Theorem 12.11 we obtain that every solution a £ V of 
P[Di,V 2 ]u = f has the form 

1 1. , 1 

U = -Ml + -(Pi + l)U 2 - ~U 3 

where U\, u 2 and u 3 satisfy the equations P 2 [T>i, T> 2 ] P 3 \P\, T> 2 ] U\ = f, 
P x \V u V 2 \P z \D u V 2 \u 2 = f and P x p x , V 2 ) P 2 [V X , V 2 ] u 3 = f. Note 
that such expression for u is not generally unique. 

Summarising, we have reduced the original equation P[Vi, V 2 ]u = /, 
see ( ITOl) . for V 1) T> 2 e End(V) to the system of four equations 

P^DuV^Pjpi^Uk = f where {i,j,k} = {1,2,3}, % < j, 

{ ) P Q [V 1 ,V 2 ]u = f. 

If we put T)\ := ^ and T> 2 := the original problem has order 6 and 
the resulting system (TT2l has order 3. 

4. Systems of polynomial equations 

The notion algebraic decompositions from [6] summarised in Section 
[2] can be applied also to systems of equations of the form (PQ) with 
commuting Pi. Here we describe one possible type of such a system to 
demonstrate power of this machinery. We will consider only the (true) 
decompositions. 
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Let us consider a /c-tuple of commuting linear operators pM £ End V 
and corresponding equations 

(13) P^u = f, [P®,P®] = 0, f E V, 1 < z, j < fc. 

The necessary (i.e. integrability) condition for existence of a solution u 
is obviously 

(14) P (i) / j = P ij) f, 1 < i < j < k. 

If, for some i, PW j s of the form jlj) satisfying (J2J), we can replace 
PW« = /W with several simpler equations using Theorem 12.11 But 
even if this is not the case, one can obtain a decomposition using an 
algebraic relation between the PWs. 

Let us consider the special case where just one equation from ( TTBl ) is 
of the form (pQ) and we do not decompose the remaining ones, i.e. we 
have the system 

Pu = PoA . . . Pei = f, [P h Pj] = 0, 1 < i, j < £, 

{ ~ l " ] R U) u = g j } [P u R ij) ] = 0, 1 < j < k, 1 < % < I 

where Pj and Rj satisfy the identity 

id v = Q P° + ... + Q e P e + S X P (1) + . . . + S k R {k) , 

(16) [Qi,Pj] = [S p ,R q ] = 0, l<i,j<e, l< P ,q<k, 
[Qi, P (p) ] = [S p , Pi\ = 0, 1 < i < i, 1 < p < k 

where P % := H{^j =0 Pi,i — 0, • • • ,£. (Note in (fl~5l l we require not only 
commutativity of the left hand sides of the equations in the systems as 
in (fT3l) but also commutativity of the factors Pj and the left hand sides 
R (j) .) The condition flU]) then becomes 

(17) R ij) f = Pg j , R {j) g i = R {i) g j for all l<i,j< k. 

Proposition 4.1 (Dual decomposition). Let us consider the system 
( f75p . Then ( fT^j is equivalent to 



(18) id v = QijPi + Qj,i p j + QljR {1) + ■■■ + QtjR (k) 

where Qi.j,Q p Lj E End(V) satisfy [Qij,P s ] = [Q 4J P (p) ] = [Q{ 3 ,R (q) ] = 
[Q P ij, P s ] = where < i, j, s < £, i ^ j and 1 < p, q < k. 

Proof. This is just a straightforward modification of the proof of Propo- 
sition EH □ 

In this setting, we obtain an analogue of Theorem 12. II for the special 
case of decompositions. In this Theorem, we replaced the operator P 
given by (pp) satisfying (j2J) with the system ((31) of £+1 simpler equations. 
Here we replace the system (fT5l) with a "system of simpler systems" as 
follows. 
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Theorem 4.2. Let V be a vector space over a field ¥ and consider 
P, : V -> V as in ( TPjjj ifte factorisation giving the decomposi- 
tion (Gi|). //ere and &e/ow we assume the range 1 < j < k. Let us fix 
f,gi G V. There is a 1-1 relationship between solutions u E V of [T5\) 
and solutions (no, • • • , Ut) E © £+1 V o/ £ne problem 

P u = /, P (1) n = PV , ■ • • , P (fc) n = P°« fc 

(i9) ; 

Pen = f, P ( % = P l g\ • • ■ , P (fc V = Py • 

Writing Vp g for the solution space of fTW and (for i = 0, • • • ,£) 
V/' s for the solution space of the system corresponding to the ith line 

in ( T71J . where g = (a 1 , . . . , g k ) T . The map F : Vp g — > X- =0 V/' g «s 
given 6?/ 

n ^ (P°n, ■ ■ • , P £ n) , 
with inverse B : X i=0 V/' g — > Vp given by 

i=l j=k 

(n , • • • , ut) h-» + • 

i=0 j=l 

On V we have B o P = id v f, e , while on the affine space X i=0 V/' g we 
/tawe P o P = id ^ t . „ . 

X i=0 V / 

Proof. Suppose it is a solution of (fT5lh Then PjP J n = Pn = / and 
also R^P' l u = P\R^u) = Py. Hence Pn is a solution of $13). For 
the converse suppose that (no, ■ • • , ue) is a solution of (1T91 ) and write 

M : = EES Qt«i + Ej=f s^'- Then 

i=£ j=k i=l j=k 

i=0 j=l i=0 i=l 

where we have used Pg- 7 = By>f from (fTTj) and then (TTBl) . Further 

«=f p=fe 



#c?) u = ^pWg iU ^ + + ( Y R U) s P g p ] 

i=0 j¥=P=± 
i=l p=k 

= (^Q.pV) + P (i) % j + ( Y S P P (j V) = g j 

i=0 3t^P=1 

where we have used R^m = P l g j from (pj) and R^g p = R^g j from 
(TT71) in the middle equality and (TTBl) in the last one. 
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It remains to show that F and B are inverses. Clearly u G Vp g 
satisfies 

i=l j=k i=£ j=k 

u = u 

i=0 j=l i=0 j=l 

since gi = Ry'u according to (TT5T) . Thus we obtain B o F = id v f, s . To 
compute the opposite direction we need the r th component of FB(u , . . . , ut) 
for (uq, . . . , Ut) G x| =0 V/' g . This is equal to 

i=l j=k i=l j=k 

P r Y,Qi u i + pr Y. S i9 j = S Q i P r u i +Q r F r u r +^2s j R^Ur = u r . 

i=0 j=l r=£i=0 j=l 

since P r Ui = P l u r foir^i (the consistency condition given by PjMj = 
/ for % = 0, ...,£) and P r g j = R^u r from (EE9]). Hence F o B = 

id Y e v /,g. □ 



5. Higher symmetries of operators 

For a vector space V and a linear operator P : V — » V, let is say that 
a linear operator 5 1 : V — > V is a formal symmetry of P if PS = S'P, for 
some other linear operator S' : V — » V. Note that S 1 : ftf(P) — > J\f(P). 
In [6] we called operators with the latter property "weak symmetries" 
and discussed the structure of the algebra of these in relation to symme- 
tries and related maps for the component operators Pi. We show here 
that although formal symmetries are defined rather differently similar 
results hold using our general tools as discussed above and in [6j. For 
the case of P a differential operator the formal symmetries agree with 
the "higher symmetries" considered in [3] and we thank Mike Eastwood 
for asking whether the ideas from [6] might be adapted to deal directly 
with what we are here calling formal symmetries. 

Consider the case of an operator P = PqP\ ■ ■ ■ Pe with a decomposi- 
tion 

(20) id v = Q P° + ■■■ + Q e P e , 

i.e. (J2J) with a := {J C L \ \J\ = 1}. Then as commented in Remark 
12.21 upon restriction to J\f(P), the operators 

Pr.^QiP', 1 = 0,1,...,* 

are projections onto Af(Pi). This was the critical object used to discuss 
weak symmetries and their decompositions in [6j. Here we see that it 
plays a similar for formal symmetries. 

Now note that if S is a formal symmetry of P then PriSPri is a for- 
mal symmetry of Pi. More generally, using the assumed commutativity 
as in ([2D, we have 

Pi(Pr t SP rj ) = Q i (PS)Pr j = Pr l (S'P)Pr J = iO ; S' Pr ,F') I) 
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so Sij := PviSPrj linearly maps Af(Pj) -> A/"(P). (In fact P*£ would 
suffice (see the remark below), we use PriSPrj for the link with [6].) 
But we may view the property -PjSV, = S'^Pj (with S'^ any linear en- 
domorphism of V) as a generalisation of the idea of a formal symme- 
try. If we have such a generalised formal symmetry for all pairs 
i,j G {0, 1, • • • ,£} then note that for each pair i, j we have 

PS tJ Pr 3 = I-PS,,!*,-, = P % S' l3 P 3 Pr 3 = (P'S'^P; 

SijPrj (and hence also PriSyPrj) is a formal symmetry of P. Thus 
the decomposition of the identity (T20l) allows us to understand formal 
symmetries of P in terms of the generalised formal symmetries of the 
component operators Pj, i — 0, 1, • • • , I. 

Remark 5.1. This result for formal symmetries follows the Theorem 
4.1 in [6] where weak symmetries are treated. The decomposition of 
the identity (1201 ) plays the crucial role in this theorem. Since, upon 
restriction to Af(P), the Pr^ are projections, the formulae above have 
a straightforward conceptual interpretation. However there is, in fact, 
an even simpler relationship between formal symmetries S of P and 
generalised formal symmetries S^, i,j G {0, ...,£}. We simply put 
:= P l S\j l /(p.) for a formal symmetry S and S := S^-P- 7 for a gener- 
alised formal symmetry S^. 

Using a factorisation from [5], this observation enables a treatment 
of the higher symmetries of the e.g. the conformal Laplacian operators 
of [7] on conformally Einstein manifolds. In particular an alternative 
approach to the higher symmetries of the Paneitz operator which is 
alternative to that in [4]. (In fact in [4] they consider only the square 
of the Laplacian on Euclidean space but by conformal invariance this 
may alternatively treated via the Paneitz operator on the sphere.) This 
will be taken up elsewhere. 
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